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Quantum light in the turbulent atmosphere 
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Nonclassical properties of light propagating through the turbulent atmosphere are studied. We 
demonstrate by numerical simulation that the probability distribution of the transmission coefficient, 
which characterizes the effects of the atmosphere on the quantum state of light, can be reconstructed 
by homodyne detection. Nonclassical photon-statistics and, more generally, nonclassical Glauber- 
Sudarshan functions appear to be more robust against turbulence for weak light fields rather than 
for bright ones. 
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Nonclassical properties of quantum light have been 
of great interest from the viewpoint of fundamentals of 
quantum physics and for a variety of applications, such 
as quantum information processing and quantum metrol- 
ogy. Special knowledge on the propagation of quan- 
tum light through the turbulent atmosphere is required 
in the context of implementations of quantum cryptog- 
raphy for communication channels between earth-based 
stations Jj and between satellites and Earth-based sta- 
tions 0. The theory is well established for the prop- 
agation of classical light through the atmosphere, see, 
e.g., 3, including phenomena such as beam wander, 
beam spreading, scintillations, degradation of spatial 
coherence, and others. However, nonclassical proper- 
ties, such as sub-Poissonian statistics of photocounts [J, 
quadrature squeezing 51, non-positivity of the Glauber- 
Sudarshan P function [6., JJ , and entanglement [8.] have 
been little studied in the context of the propagation of 
light through the turbulent atmosphere. 

Due to the occurrence of random fluctuations of the 
refractive index, the quantum state of light after trans- 
mission through the atmosphere cannot be presented 
by a single-mode density operator, in terms of neither 
monochromatic nor nonmonochromatic modes. The pho- 
tocounting statistics of a combination of scattered modes 
has been studied by a random modulation of the inten- 
sity [§, • This model has been further improved [III ■ 
The technique of the photon wave-function allows one 
to consider special cases of single-photon and two- 
photon [l3| states. Another approach, which describes 
single-photon states, is presented in Ref. p^ . 

In the present contribution we deal with the effects 
of an atmospheric transmission channel on any quantum 
state of light. General expressions for the quantum state 
after transmission are derived. Based on balanced homo- 
dyne detection, one may reconstruct the statistical dis- 
tribution of the transmission coefhcient through the at- 



mosphere. By repeated reconstruction of the statistical 
distribution, one can significantly reduce the atmospheric 
noise effects on the quantum state of light. 

For dealing with continuous- variable quantum states, 
we are considering the experimental setup in Fig. [1] 
The light from a source is transmitted through the at- 
mosphere and collected by a telescope (or some other 
device). Subsequently, a balanced homodyne detection 
setup is used to filter out the desired nonmonochromatic 
mode from other modes and background radiation b y an 
appropriate local oscillator, for details see [H, [H, 
The remote local oscillator can be synchronized with the 
source field by, e.g., the technique of the optical frequency 
comb [l^ . In the limit of a strong local oscillator, the dif- 
ference of photocurrents in the detectors is proportional 
to the field quadrature of the nonmonochromatic output 
mode defined by the local-oscillator pulse. Knowledge of 
the quadrature distributions enables one to get the com- 
plete information about the quantum state of the consid- 
ered output mode. For example, the method of optical 
homodyne tomography enables one to reconstruct the 
Wigner function, the photon-number distribution, mo- 
ments of the radiation field, and the density operator in 
an arbitrary representation, for a review see [iq . 
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FIG. 1: Homodyne detection of quantum light generated by 
the source, S, and propagated through the turbulent atmo- 
sphere. BS is a 50:50 beam-splitter, Di and D2 are detectors, 
LO is the local oscillator. 

We assume that the propagation of quantum light 
through the turbulent atmosphere is a linear attenuating 
process. The corresponding losses are caused by light ab- 
sorption and scattering as well as by the mode-mismatch 
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due to different shapes of the output and local-oscillator 
pulses. The Glauber-Sudarshan P function [19|] of an 
attenuated nonmonochromatic output mode, Pt (a) , is 
related to the P function of the input mode, Pin (a), as 
(cf., e.g., 0) 



(1) 



where T is the complex transmission coefhcient, with 
\T\^ < 1. An important difference between turbulent me- 
dia and other lossy systems is that the transmission coef- 
ficient r is a random variable with fluctuating phase and 
magnitude. This means that the P function of the out- 
put mode. Pout (ck), is obtained through averaging Pt (a) 
with the probability distribution of the transmission co- 
efficient (PDTC), r{T), as 



Pont (a) 
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The integration is performed over the circular area, 
|r|^ < 1. This represents the quantum-state input- 
output relation for light propagated through the turbu- 
lent atmosphere. 

The explicit form of the PDTC should be obtained 
from a theory, which considers turbulence properties of 
the atmosphere as well as specific conditions of the ex- 
periment. Since this is a complex problem, we may only 
consider a simple model for the case of small fluctuations. 
The corresponding PDTC can be obtained similar to the 
probability distribution of the intensity modulation in 
Ref. [9]. For this purpose, we consider a discrete set of 
turbulent eddies, each of them is characterized by a ran- 
dom transmission coefficient Tk- The total transmission 
coefficient is T — Y[Tk- The central limit theorem im- 

k 

plies that the PDTC is a two-dimensional distribution, 
which is log- normal with respect to the magnitude i = |r| 
and normal with respect to the phase (p — argT, 
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(3) 



Here, and ag are the mean value and the variance, 
respectively, of 9 = —Int; acp is the variance of ip; s is 
the correlation coefficient between 9 and ip. Without loss 
of generality we suppose that the mean value of (p is zero. 
This form of the PDTC can be used only for ag <^ 9 and 
a^p <C 27r. Contrary to the approach based on the random 
modulation of the intensity j^, [l^, [HI , we restrict the t- 
integration to the range < t < 1. Accordingly, the 
(^-integration is restricted to — tt < ip < tt. 

Of course, the given model [Eq. ^] will not properly 
describe the turbulence properties of the atmosphere un- 
der general conditions. Due to the lack of a general 
model, it is important to develop a method for the ex- 
perimental determination of the PDTC. Let the input 



field be prepared in a coherent state I7). In this case the 
PDTC can be expressed in terms of the characteristic 
function $out (/?), of the P function of the output state 
as 
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where and Ti are the real and imaginary parts of the 
transmission coefficient, respectively. In an optical homo- 
dyning experiment ^out can be estimated from a sam- 
ple of N photocounting difference events An^, cf., [tI-IT^. 
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The amplitude and the phase of the local oscillator are 
fixed to be r and (f — arg/3), respectively. 

To demonstrate the practical usefulness of the recon- 
struction method of the PDTC, we have performed the 
following simulation. We start with model distribu- 
tion ([21). For simplicity, it is approximated here by a 
normal distribution in the variables Tr and Ti. In prac- 
tice, observed data shall be used, so that this assumption 
does not restrict the applicability of our method. We de- 
rive the P function of the output field from Eq. ^ and 
calculate the photocount-difference distribution by using 
the corresponding integral transformation, see Ref. [2ll |. 
Now we can simulate the measured data and reconstruct 
by Eqs. (g]) and (H]) the PDTC. The result of this proce- 
dure is shown in Fig. [2l which is in reasonable agreement 
with the initially chosen PDTC even for a small sample 
of data. In real experiments, this method yields insight 
into the true statistics of the turbulent atmosphere. 




FIG. 2: (Color online) The PDTC is shown for a normal- 
distribution approximation of Eq. ((3} with 8 — 0.9, ae = 
0.2, = 0.2, s = 0.01, Ti = (a), 0.1 (b). The solid and 
the dashed lines are the initially chosen and the reconstructed 
(from only 5 x 10^ sampling events) PDTC, respectively. 

Alternatively, the PDTC can also be characterized by 
its statistical moments. Consider the matrix of normal- 
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ordered moments of the photon annihilation (creation) 
operator a [a)) for the source field, 



M„„ = Tr [pat"a"'] = J a P (a) a*"" a"" , (6) 



which completely characterizes the quantum state of a 
light mode, where p is the density operator. The mo- 
ments can be measured by balanced homodyne detec- 
tion Utihzing the input-output relation [Eq. Q], 
one gets the corresponding relation between the normal- 
ordered moments of the input and output modes. 



where 
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are the moments of the PDTC. From Eq. ([7]), the mo- 
ments of the PDTC, 
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are obtained by measuring the moments of the input and 
output fields. For example, one can use the input mode 
in a coherent state I7) such that M^^ = 7*"7™. The 
obtained moments of the PDTC allow one to determine 
the nonclassical properties of the output mode, once the 
corresponding properties of the input mode are known. 

It is known |3| that due to the atmospheric winds nar- 
row light beams are randomly deflected and wide beams 
scintillate within a certain time Tatm- For data accu- 
mulation times Tdata ^ Tatm wc cxpcct large fluctua- 
tions of the transmission coefficient T. Otherwise, for 
Tdata ^ T"atm the transmission coefficient is not signifi- 
cantly fluctuating during the time Tdata- However, for 
this scenario, V (T) is randomly changed between differ- 
ent series of measurements, separated by time intervals 
At > Tatm- The PDTC and its moments can thus be 
permanently monitored by using the method proposed 
above, which works with a small sample of data that can 
be recorded within short Tdata intervals. In this way one 
can suppress the influence of long-term atmospheric fluc- 
tuations on the quantum state of the transmitted light. 

Let us consider the transmission of sub-Poissonian 
light 3 through the turbulent atmosphere. Using 
Eq. ([7]), the Mandel parameter (cf., e.g., [20|) of the out- 
put field, Qout, can be related to the Mandel parameter 
of the input field, Qin, as 
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where r] = T*T is the efficiency. The first term resembles 
the behavior for standard attenuation. The second term 



is caused by fluctuations of the efficiency rj due to the 
atmospheric turbulence. It is proportional to the mean 
photon-number of the input field, M{" = {n)^^^. For states 
of the input field whose mean photon-number fulfills 
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the photocounts of the output mode are always super- 
Poissonian. Hence the nonclassical photon statistics of 
bright quantum light is destroyed by fluctuations of the 
magnitude t of the transmission coefficient, but it is not 
affected by phase noise. 

In other cases, the nonclassical properties of bright 
light can also become sensitive to both phase and mag- 
nitude noise. In the most general case a given quantum 
state is nonclassical if its P function is not positive def- 
inite 0. For weak turbulence, the PDTC has a strong 
maximum at T = Tq. Hence the input-output relation 
[Eq. ^] reads in the first-order Laplace approximation 
as 
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(12) 

where 7 = (a)in is the displacement parameter of the 
input field and the PDTC is taken in the Gaussian ap- 
proximation [2^ . If the minimum eigenvalue of the co- 
variance matrix of the scaled PDTC in Eq. obeys 
Amin > 2, Pout (a) represents the Husimi-Kano Q func- 
tion [2^ of the displaced input field combined with a 
Gaussian noise, which is always non-negative. Based on 
the above approximations, we derive that for any input 
state with 



|7| > 2e' 




the corresponding output state is classical. For simplicity 
we have considered only real 7 and = 0, the general- 
ization to complex 7 and arbitrary (p is straightforward. 

As an example we consider the P function of displaced 
single-photon-added thermal states (SPATSs) Q, 



Pin (a) 
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(14) 

where nth a-nd 7 are the mean number of thermal photons 
and the coherent displacement amplitude, respectively. 
We compare the P functions of the displaced SPATS 
for two cases: (i) the standard attenuation with a fixed 
transmission coefficient T in Eq. ([T]), and (ii) transmis- 
sion through the turbulent atmosphere as described by 
model ([3]). The mean values of the transmission coeffi- 
cient are chosen to be equal in both cases. The standard 
attenuation does not destroy negativities of the P func- 
tion. For small displacements, the situations in the cases 
(i) and (ii) remain similar to each other. However, with 
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increasing displacements 7 the atmospheric turbulence 
destroys the nonclassical effects, even when they survive 
for standard attenuation, see Fig. [3l 
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FIG. 3: (Color online) P function of the displaced SPATS, for 
Im(a) = 0, fith = 1.11, 7=1 (a), 7 (b), 20 (c). The dashed 
lines show the standard attenuation for T = e"'''^ « 0.7408. 
The solid lines represent the model ^ for 6 — 0.3, ae — 0.1, 
= 0.14, and s = 0.01. 



In conclusion, we have studied the effects of atmo- 
spheric turbulence on the quantum properties of light. 
It has been shown that the probability distribution of 
the transmission coefficient can be experimentally recon- 
structed by homodyne measurements. Based on such a 
method, one may predict the general turbulence effects 
on any quantum state of light. By repeated short-time 
monitoring of the PDTC, one can suppress the disturb- 
ing effects of long-term fluctuations on the quantum state 
of light. Balanced homodyne detection also allows one 
to reduce the effects of background radiation, which is 
useful for quantum communications under day-light con- 
ditions. We have shown that the nonclassical effects of 
bright light fields can be more fragile against turbulence 
than for the case of weak fields. A nonclassical photon- 
statistics is destroyed by turbulence if the mean photon 
number exceeds a critical value. General nonclassicality 
of the P function in the output channel is sensitive to the 
mean coherent amplitude of the input radiation. 
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